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ABSTRACT 


This thesis provides closed form expressions for the angular distribution in 
azimuth and elevation planes for a geometrically based single bounce spheroid model. 
The geometry of the spheroid is defined by the semi-major axis a and the semi-minor 
axis b. The other parameter of interest in the model is the distance D between the base 
station and the mobile station. The latter is assumed to be at the center of the spheroid. 
The mobile station is assumed to be the transmitter, while the base station is the receiver. 
This thesis investigates the effects of the above parameters on the angular distribution of 
the received waves. Important parameters such as the r.m.s angle spread in azimuth and 
elevation plane are calculated from the p.d.f. expressions derived. The behaviour of these 
r.m.s angle spreads versus the ratio a/D or b/D respectively is also investigated. 
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EXECUTIVE SUMMARY 


This thesis provides closed forms which describe the statistics of the angular 
distribution in azimuth and elevation planes for geometrically based single bounce spheroid 
model. This spatial model can be generated by revolving an ellipse, which lies on the xz 
plane of a 3D Cartesian coordinate system, about the x-axis. This ellipse has semi-major 
axis a and semi-minor axis b, respectively. The base and the mobile stations are located in 
the xy plane with the mobile station being at the center of the spheroid and the base station 
at the origin of the coordinate system. It is assumed that the mobile station is transmitting 
and the base station is receiving. The base and mobile stations are separated by a distance 
D. 

The scattering spheroid model is applicable to a macrocell environment in which the 
received multipath signals at the base station are originated after a scattering from the 
surrounding environment about the mobile station. Assuming that the scatterers are 
uniformly distributed with a constant scatterer density function in the volume V as defined 
by the geometry of the spheroid, we derive the marginal angle of arrival (AOA) probability 
density functions (p.d.f). in the azimuth and elevation plane. 

Using the derived closed forms for both p.d.f.s we generate plots to examine how 
the parameters a, b and D affect the angular distribution of the received waves. From the 
marginal AOA p.d.f. in the azimuth plane we conclude that: 
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1. 


The higher the ratio a/D, the lower the probability for the received 
multipath signals at the base station to be confined to a small angular 
region centered about the axis which connects the base and the mobile 
stations. 

2. The angular distribution in the azimuth plane is independent of the length 
of the vertical dimension b of the spheroid. 

3. Comparing the results of the spheroid model for small values of b with 
those of the circular scattering model presented in [Ref. 3], it was 
observed that there exists a significant difference between them. 

Additionally, the effect of the a/D ratio on the r.m.s angle spread in azimuth plane is 
discussed. From the plot results of the r.m.s angle spread in azimuth plane versus a/D, 
we conclude that for values of a/D up to 0.5, there exists a linear relation. However, for 
higher values of a/D, the relationship changes to a nonlinear curve with a positive 
gradient. 

Regarding the marginal AOA p.d.f. in the elevation plane, the following 
conclusions were drawn from the plots: 

1. Similar to the p.d.f in the azimuth plane, the higher the ratio b/D the lower 
the probability for the received multipath components at the base station to 
be restricted in a small angular spread around the axis, which connects the 
base and the mobile stations. 


2. 


The r.m.s angle spread which also depends on the semi major axis a, has 
a linear dependence with respect to the b/D ratio. 
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I. INTRODUCTION 


A. BACKGROUND 

The outdoor radio environment contains a large number of obstacles, such as 
buildings, cars, hills, etc, which obstruct, reflect and scatter the wireless radio signals. As 
a result, the transmitted signal arrives at the mobile station or at the base station from 
many different angles and times of arrival; this is known as multipath. Another common 
phenomenon in cellular systems is the co-channel interference, and it is created because 
in a certain geographic region(cluster) there are other subregions(cells) that utilize the 
same set of frequencies. Multipath fading and co-channel interference can increase 
significantly the bit error rate and they don’t allow high bit rate data services in wireless 
communications systems. The presence of multipath and the motion of the receiver or the 
transmitter causes the received signal envelope to be time-varying and experience a 
Doppler shift that depends upon the angle of arrival of the incoming waves. Concepts 
such as angular spread, time delay spread, and Doppler spread have to be taken in 
consideration during a wireless communication system evaluation. In particular, they 
should be taken into consideration in the design of smart antennas used for diversity, 
beamforming and emitter localization applications. 

The multipath effect alone can be countered by applying antenna spatial diversity 
which implements an antenna array consisting of multi-elements. Spatial diversity is an 
arrangement whereby signals from multiple elements are combined to result in a signal 
that fades less rapidly than the individual element signals. Ideally this can be 
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accomplished if received signals at the two elements are decorrelated. The degree of the 
decorrelation between the two received signals by the two elements can be expressed by 
the spatial cross-correlation coefficient which is a function of the joint p.d.f. L e (</>,0) 

describing the angular distribution of the incoming waves in the three-dimensional space 
at the receiver. 

The interference limited factor can not be overcome by using spatial diversity. 
Switched beam arrays can be used with some success if the angle of arrival of the 
interfering signal is not near the angle of arrival of the desired signal [Ref. 1], Only the 
beamforming arrays have the potential to completely suppress the co-channel interference 
(if the number of elements exceeds the number of co-channel interferes) improving the 
wireless communication system performance. The performance of these type of smart 
antennas depends on the spatial characteristics of the wireless radio channel. The term 
“spatial characteristics of the radio channel” is used to mean the statistics of the received 
signal envelope in the spatial domain. These spatial characteristics may not always be 
available from measurements. In such a case one needs to devise geometric models to 
depict the scattering process in the spatial domain. Another area where knowledge of 
angular distribution is valuable is in the determination of the power spectral density [Ref. 
1]. Therefore the angular distribution of the received multipath components is important 
in evaluating the performance of radiowave link between the transmitter and the receiver 
in wireless communications. 
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Spatial channel models have been created and applied in order to be able to derive 
the statistics of the received signal envelope for both macrocell and microcell 
environments. The term macrocell enviroment assumes that the scatterers surrounding 
the mobile station have the same height or are higher than the antenna of the mobile. 
This results in a received signal at the mobile station which arrives from all directions 
after bouncing from the surrounding scatterers. In contrast, the base station’s antenna is 
located higher than the surrounding scatterers. Therefore the received multipath signals 
at the base station originate after scattering from the area near to the mobile station 
[Ref.2], In microcell enviroment both base and mobile station are deployed in the same 
height as the surrounding scatterers. This implies that the scattering process takes place 
near to the base and mobile station. The characteristic difference between the types of 
environments is that the angular spread of the received signals at the base station is 
confined to a small region for macrocell environments. In the microcell enviroment this 
region tends to be larger. 

In Table 1, which is a modified version of a table from [Ref 2.], are listed some 
representative spatial channel models. 
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SPATIAL CHANNEL 

MODEL 

DESCRIPTION OF THE 
MODEL 

APPLICATIONS 

1. Lee’s Model 

Scatterers are evenly spaced 
on a circular ring about the 
mobile at a distance R. Each 
of the scatterrers represents 
the effect of many other 
scatterers and hence is 
referred to as an effective 
scatterer. The base station is 
deployed a distance D from 
the mobile station. 

Predicts correlation 

coefficient using a discrete 
angle of arrival (AOA) 
model. A modified version 
of this model can account 
for a Doppler shift, 
assuming that the scatterers 
are moving with an angular 
velocity. 

2. Discrete Uniform 

Distribution 

A model similar to Lee’s 
model with the N scatterers 
evenly spaced over an 
AOA range. 

Predicts correlation 

coefficient using a discrete 
AOA model. The 

correlation predicted by this 
model falls off more quickly 
than the correlation in Lee’s 
model. 

3. Geometrically Based 
Circular Model 

In this model it is assumed 
that the scatterers are 
uniformly deployed within a 
radius R about the mobile. 
The locations of the 
scatterers are defined by a 
spatial scatterer density 
function. 

It is suitable for macrocell 
environments and predicts 
the joint time of 

anival(TOA) and AOA and 
marginal density functions 
at both base station and 
mobile station. 

4. Geommetrically Based 
Elliptical Model 

Scatterers are uniformly 
distributed in an ellipse 
where the base station and 
mobile station are the focii 
of the ellipse. 

It can be determined the 
marginal AOA, TOA, the 
joint TOA and AOA density 
function. It is applied in 
microcell enviroment. 

5. Gaussian Wide Sense 
Stationery Uncorrellated 

Scattering (GWSSUS) 

N scatterers are grouped 
into clusters in space such 
that the delay differences 
within each cluster are not 
resolvable within the 

transmission signal BW. 

Provides an analytical 
model for the array 
covariance matrix. 
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6. Gaussian Angle of Special case of the Provides an analytical 
Arrival (GAA) GWSSUS model with a model for the array 

single cluster and angle of covariance matrix, 
arrival statistics assumed to 
be Gaussian distributed 

_ about some nominal angle. _ 

Table 1. Description of some representative spatial channel models. “From Ref. [2].” 

Figure 1 through 4 demonstrate the first four spatial channel models, listed in 
Table 1. 



Figure 1. Lee’s Model “From Ref. [2].” 






Figure 2. Discrete Uniform Distribution Model “From Ref. [2]”. 



Figure 3. Geommetrically Based Circular Model “From Ref.[2]”. 



















Figure 4. Geommetrically Based Elliptical Model From Ref.[2]”. 


Every spatial channel listed in Table 1 has one common feature; the received 
signal at the base station is assumed to arrive from a two dimensional plane connecting 
the tips of the transmitter and the receiver. Hence the angle of arrival (AOA) includes 
only the azimuth angle information. 

In this thesis we consider a 3D spheroid model as shown in Figure 5 This 
spheroid can be generated by revolving an ellipse, which lies on xz plane, with semi¬ 
major axis and semi-minor axis a and b, respectively, about the x-axis. The base station 
and the mobile station are located in the xy plane with the mobile station at the center of 
the spheroid, and the base station at the origin. It is assumed that the mobile station is 
transmitting and the base station is receiving. The distance between the two terminals is 
D with the requirement that the semi-major axis of the spheroid be less than or equal to 
















B. OBJECTIVE 

Assuming that the scatterers are uniformly distributed around the mobile station 
as defined by the geometry of the spheroid, the objective is to derive the marginal 
probability density function (p.d.f.) angles of arrival (AOA) in the azimuth and in the 
elevation planes at the base station. After extracting these quantities, we investigate how 
the angular spread in azimuth and elevation plane are related with the dimensions of the 
spheroid and the distance between the transmitter and the receiver, as well as how they 
are correlated with each other. 
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II. PROBLEM FORMULATION 


A. INTRODUCTION 

In this chapter we consider a single bounce scattering spheroid model. The model 
can be used in a macrocell environment where the base station antenna height is 
relatively high. As a result there will be no scattering from scatterers near the base 
station. Figure 6, shows the spheroid scatterer density geometry and the notation used to 
derive the marginal p.d.f.s. As is shown in the Figure 6, the location of each scatterer is 
defined by the Cartesian coordinates (x,y,z) or the spherical coordinates (r b ,<p b ,0 b ), or 

{ r b>0b’fi)’ where /3 represents the complementary angle of 6 b in elevation plane. 
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Figure 6. The spheroid scatterer density geometry. 

The method for the derivation of the joint TOA/AOA p.d.f. and marginal AOA p.d.f. for 
both azimuth and elevation planes is based on the approach presented in [Ref. 3]. In the 
Figure 6, it is assumed that scatterers are distributed randomly in the volume V of the 
spheroid according to the spatial scatterer density function f xyz (x,y,z ). It is further 

assumed that the waves arrive directly at the base station after undergoing scattering only 
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once from within the volume V. Multiple scattering of waves is not considered in this 
thesis. 

Before proceeding with the determination of the p.d.f.s in the azimuth and the 
elevation plane we would like to present the formulation for the joint TOA/AOA p.d.f 
and marginal AOA p.d.f in the general case. In section B we start with the derivation of 
the joint TOA/AOA p.d.f.. In section C, we show formulation for the marginal AOA 
p.d.f.s. In the following sections D through G we derive the marginal AOA p.d.f in 
azimuth and elevation plane for the spheroidal model. 

B. DERIVATION OF THE JOINT TOA/AOA P.D.F.’S 

Figure 7, shows a plane wave impinging on a base station after scattering from an 
obstacle, represented by a point ‘S’ in space with the ordered rectangular triples (x,y,z). 
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The scatterer probability density function f x>y>z (x, y, z) can be also expressed 

terms of spherical coordinates systems (r b ,0 b ,6 b ) as follows. The equations relating the 
spherical coordinates to the Cartesian coordinates are: 


r b 

o b 


-V 


2 2 2 
x z + +z 


■ cos 


-ii 


^x+y+z 2 


( 1 ) 

( 2 ) 
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fa = tan : (j) 

(3) 

x = r b sin 0 b cos fa 

(4) 

y = r b sin 6 b sin fa 

(5) 

z = r b cos0 b 

(6) 


where (x, y,z) and (r b , fa, d b ) denote the location of the scatterer in Cartesian and 
spherical coordinates respectively. The joint p.d.f f h ^ ^ ( r b’fa>@b) can be found using 
the well known transformation [Ref. 4]: 


x = r b sin 0 b cos <j> b 

fr b 4b A ( r b ’fa,Ob) = J ( r b . fa > fa)fx,y,z (•*. y>z}y = r b sin d b sin (j) b 

z = r b cosd b 

where J(r b ,fa,0 b ) is the determinant of the Jacobian transformation given by: 


(7) 


j{ r b’fa’&b) = 


dx dx dx 


dr b d<p b dO b 
dy dy dy 


dr b d(j) b d9 b 
dz dz dz 


dr b d(j> b dQ b 


= rjj sin 6 b , 0 <6 b <n 


( 8 ) 


Substituting the equation (8) into (7) we get for the scatterer probability density function 
fr bAA i r b’fa’@b) w ' tb respect to spherical coordinates: 


fr b 4bA ( r b’fa’ 0 b) = r b sin fafx,yA r b sin0 b cos fa,r b sin 6^ sin fa,r b cos0 b ). (9) 

If we want to express the scatterer probability density function f h ^ e ( r b’fa’&b ) 
with respect to (T, fa,6 b ), where T denotes the time delay of the multipath component, 
we have first to find the relationship between r b and the T. Applying the law of cosines 
to the triangle BS,S,MS we get the following equation: 
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r 2 = r 2 + D 2 -2r b DcosS 


( 10 ) 


where the cos <5 is called the ‘cosine direction’ given by: 

cos S = sin 6 b cos (j> b . ( 11 ) 

The delay corresponding to the total path traveled by the wave is: 

T = = 7 [ r b + V? +D 2 + 2r b D sin 0 b cos <p b ) (12) 

where ‘c’ is the speed of propagation of waves. 

Squaring both sides of the equation ( 12 ) and solving for gives: 

(13) 

The joint TOA/AOA p.d.f. f z ^ e ^{r,^ b ,0 b ) can be found using the Jacobian 
transformation [Ref. 4] according to: 




o 2 -A 2 


y i i _ 

^ 2(Z)sin0 6 cos ^,—rc) 


(14) 


where is the determinant of the Jacobian tranformation given by: 




2(£>sing fc cos^, -re) 2 
0 2 c+r 2 c 3 -2-rc 3 £>sin ^ cos( 


(15) 


Substituting the equation (15) into (14) we get for the joint TOA/AOA p.d.f. 


A** ( 16 ) 


Using the last equation (16) and the equations (9) and (13) we can express the joint 
TOA/AOA p.d.f. f T 0 bdb {t,0 b ,6 b ) in terms of original scatterer probability density 


function f xy z (x, y,z) as follows: 


16 



/rA,&( f 4.^) = sin^ 


(P 2 -r 2 c 2 ) 2 (Z) 2 c+r 2 c 3 -2g 2 Psin^cos<i 
8(/>sin 6 b cos (fy, -nr) 4 


- Sin <9* cos <j> b , r b sin d b sin <f > b , r b cosd b ) • 


(17) 

Note that the equation (17) expresses the joint TOA/AOA p.d.f. f z ^ 0 {z,(p b ,6 b ) as 
observed at the base station in terms of an arbitary scatterer probability density function 

fx,y,z( X >y> Z )- 

When the scatterers are uniformly distributed within an arbitrarily space A with a 
volume V, then the scatterer probability density function f xyz {x, y, z) is given by: 


fx,yA x *y>z) = 


y,x,y,ze A 
0 , otherwise 


(18) 


If we consider this special case, the joint TOA/AOA p.d.f. f T ^ & (t,(p b ,6 b ) can be 


written: 


fr,(p b ,e b ( r > s in 6b 


(Z) 2 -r 2 c 2 ) 2 (P 2 c+r 2 c 3 - 2a: 2 Dsin^cos4) 
8 V (D sin 6 b cos <p b -tc ) 4 


where the range for the values of T, (j) b , 6 b has to be defined. 


(19) 


C. DERIVATION OF MARGINAL AOA PROBABILITY DENSITY 
FUNCTION 

The marginal AOA p.d.f. with respect to spherical angle coordinates could be 
found by integrating the joint p.d.f. f rh ^e b {r b ,(p b ,6 b ) , as given in equation ( 9 ) with 

respect to r b over the range \{<t> b ,0 b ) to r b {<p b ,6 b ) , and then either with respect to 

sM over the range 6 b {(J) b ) t0 6 bM or with respect (j) b (d b ) over the range 

<t> bi (0 b ) to <j) b ^ ($ b ) depending on whether the marginal A.O.A p.d.f. is to be found in 
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azimuth or elevation plane respectively. The terms r b {<p b ,0 b ) and r b (<p b ,0 b ) are the 
lower and the upper limits of integration and depend on the azimuth angle </> b and the 
elevation angle 6 b . The values of these limits can be found from the intersections of a 
straight line 0 b = constant, <j> b = constant with the scatterer volume. The limits of 
integration in 6 b , ^and^for the marginal AOA p.d.f. in azimuth plane can be 
defined as the smallest and the largest 9 b -values that bound a region R, which is the 
intersection between a plane at a specific <j) b value and the scatterer region volume. The 
above statement can be expressed in mathematical expressions as follows. 

The marginal AOA p.d.f. (<p b ) in azimuth plane is given by: 

4 (& ) = C'u) J fr b &A 4 )dr b dd b = 

hK b) r^A) 

0 b 2 (<Pb) r b2 (fa’ 0 b) 

J J r b Sin 6 b f x y z (r b s\x\0 b cos<p b ,r b sir\0 b sin <p b , r b cos0 b )dr b d0 b . (20) 

% ) % {0b A) 

The marginal AOA p.d.f. f ff (d b ) in elevation plane is given by: 

4 ( 0 b ) = Jwf? J fr bA ,e b { r b’0 b ’@b )dr b d(j) b = 

^ b> r^A) 

0b 2 i e b) r b 2 ( ( PbA) 

J J r b sin 0 b f x y z (r b sin 0 b cos <p b , r b sin 0 b sin (p b ,r b cos 0 b )dr b d<j> b . (21) 

0b\ (&b ) r b[ (0b A ) 

The limits for 0 b integration, which depend on the elevation angle 0 b , can be 
determined from the smallest and the largest azimuth angles that bound the region of 
intersection between the cone 0 b = constant and the scatterer volume. 
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Equations (20) and (21) express the AOA marginal p.d.f.’s for any scatterer 
probability density function f xyz (x,y,z). When the scatterers are uniformly 

distributed in a volume with a constant probability density function f x z [x,y,z), the 
marginal AOA p.d.f. (<p b ) is given by: 

%(8 b )r bl %(<h) 

4 (*») = £ / J'i 2 sin^'i d ^=w Jr4%A)A 

°b\ i e b ) \ (h fib ) &b ( 0 *) 


( 22 ) 


and the marginal AOA p.d.f. f 0 (d b ) is given by: 


0*2 & ) r *2 ( 0 * '®b ) Vbi \&h ) 

fe b (0 b )=V / JV sin 0 b dr b d</> b J fo* (</> b ,6 b )- , 6 b )]sin d b d(l) b 

\{9b ) r b l (M) \{&b) 


(23) 


D. MARGINAL AOA P.D.F. IN AZIMUTH PLANE FOR THE SPHEROID 
SPATIAL CHANNEL MODEL 

Before we proceed to find the marginal AOA p.d.f. in the azimuth plane, we have 
to define the equation which describes the geometry of the spheroid in Cartesian and 
spherical coordinates. The spheroid in Cartesian coordinates is given by: 


a 2 b 2 


(24) 


and in spherical coordinates is given by: 

(a 2 cos 2 6 b + b 2 sin 2 6 b )r 2 - 2 Db 2 r b sin 6 b cos <j> b + b 2 (d 2 - a 2 ) = 0. 


(25) 
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Instead of the polar angle 0 b , its complementary ft — y — 0 b will be used hereafter for 
simplicity. Equation (25) can be expressed as: 

- B(fi)r b cos <j) b + C = 0 (26) 

where, A(j3) = a 2 sin 2 j3 + b 2 cos 2 P, B(fi) = 2Db 2 cos^, C = b 2 (D 2 - a 2 ) and 
A(/?) = B 2 (J3)cos 2 0 b -4A(J3)C 

= 4 b 2 [a 2 (a 2 sin 2 j3 + b 2 cos 2 /?)- D 2 a 2 sin 2 /? - D 2 b 2 cos 2 sin 2 <f> b \ 

The two roots of the above equation are: 



(27) 

In terms of the angle (3 , the marginal AO A p.d.f. (<%) in azimuth plane expressed 
as: 


fo b (<Pb) = w J Vl {<t>b >/?)- r l {0b , P )]cos pdp (28) 

*Amax {0b ) 
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z 



Figure 8. The intersected region R between a plane (j) b = constant 
and the spheroid. 

Figure 8 shows the limits of the angle ft for a specific value of (f ) b . The limits of 
integration in the elevation angle can be determined from the values of the angle /3 that 
bound the region R, which is produced by the intersection between the plane <j) b =constant 
and the spheroid. In our case the cross section of the spheroid with a plane <p b = constant 

/ 9 9 9 

is an ellipse with semi-major axis equals to yja - D sin (p b and semi-minor axis 

equals to - D 2 sin 2 (f> b and the center located at Dcos<p b . For a general angle j3 , 

the straight line from the origin intersects the region R at two points and r b ^. At the 

two boundary points the roots and r hi of the quadratic equation (26) in r b coincide. 

This happens when the discriminant of the quadratic equation is set to zero. Carrying this 
out, the upper limit /? max is given by: 
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/^max tan 


b_ la 2 -D 2 -sii 
a V D 2 -a 


•sin 2 <Pb 


(29) 


This equation is valid for (f> b in the range of sin ' (--g-)< <j> b < sin _1 (fr)- These 
correspond to the values of (j> b that bound the spheroid in the xy-plane. Taking all of 
the above into consideration, the marginal AO A p.d.f. in azimuth plane is given by: 

Am <A> 

ftS<t>b) = W Jcos p[rl(t b ,P)-rl(j b ,PW0 , sin"(-f)<4<sin _1 (f) . (30) 


The evaluation of the integral as defined in the equation (30) is straightforward. 
The following relations simplify some of the terms: 

= (% (0b’P) ~ % (<t> b , P))(r^ (4> b , p) + r h (<f) b , P)r h (<p b ,0) + r* (<p b , (3)) . (31) 


r h (<p b ,P)-r bt (<p b ,^) = -^j 


MA,P) + r b2 (<P b ,/3) = ^- 


rfat, fi) + %(<f> b , fi) = 


B 2 (J3) cos 2 


A Z (I3) 


_2C_ 

A(/3) 


r bS<Pb’0)- r b 2 (<Pb’fi)=ifyj 


Using equations (31) through (35) in equation (30) we get: 

Anax , . 

4 Wb) = w Jk + r by r h 2 + r b 2 jcos fidfi - 

o 

— _2_ f Vg 2 (yg)cos 2 ^-4A(^)C /g 2 (/?)cos 2 l>„ c 1 /U/? 

3 V J A(£) \ A 2 (/3) A(p)p OS P a P- 

0 


(32) 

(33) 

(34) 

(35) 


(36) 
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At a specific value of (f> b when /? = /? max the discriminant becomes 
Af/L*) = lens 2 0 b - 4A(fi max )C = 0. Hence, writing 

a(0=aG9)-a(O 

= [B 2 (/3)cos 2 </>„ - 4A(/3)C] - [B 2 (fi m )cos 2 </>„ - 4A{^ )C], 

we see that 


A(/?) = 4b 2 [a 2 [a 2 sin 2 fi + b 2 cos 2 /?)-Z> 2 a 2 sin 2 /3-D 2 b 2 cos 2 /?sin 2 <f> b J- 

-4fc 2 [a 2 (a 2 sin 2 /? max +& 2 cos 2 /? max )-Z> 2 a 2 sin 2 $ naJ( - D 2 b 2 cos 2 fi max sin 2 $J, 


which gives: 

A = 4b 2 [sin^ - sin 2 /?Ja 2 (l) 2 -a 2 )+b 2 (a 2 - D 2 sin 2 ^)]. 

Using this the integral in equation (36) can be written as follows: 



16£>"fr 5 cos 2 [D 2 -ar\rb 2 [a 2 -D 2 sin 2 <p b ) 
— 


■'max 

I 


cos 3 fi<Jsin 2 -sin 2 P ip 
(a 2 sin 2 f$+b~ cos 2 y?) 3 P 


4(D 2 - a 2 )t. 3 V« 2 ^ 2 -« 2 ^ 2 (a 2 -D 2 sin 2 0j 

3V 


Anax _ 2 „ . i _ 

f COS/?VSin Anax~Sm 2 ^ rg 

J (a 2 sin 2 (3+b 2 cos 2 /?) 2 


(37) 


(38) 


(39) 


The first integral in equation (39) can be simplified as: 

V cos 3 ^ysin 2 /? ma) ,-sin 2 / g _ ) _ cos^(1-sin 2 yg^sin 2 -sin 2 

' (a 2 sin 2 ^+Z> 2 cos 2 ^) 3 J (a 2 sin 2 >0+6 2 cos 2 y9) 3 ^ 


(40) 


On using the substitution t = sin /?, and r max = sin /? max 

/? ? X cos^Ct-sin 2 ^)Vsin 2 ^ max -sin 2 ^ „ _ ? T (l-r 2 X/rg Lax -f 2 
q (a 2 sin 2 p+b 2 cos 2 /?) 3 J [ & 2 + ( a 2_ £) 2^2] 3 


the last integral becomes: 


(41) 
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Hence the term l—t can be written as follows: 

Using this in equation (41) we get: 


J [»M- 2 -» 2 W 


dt = 


‘max 

= _ai_ f 

a 2 -b 2 J 


V^max \ 


N« 2 -* 2 W 


dt 


‘max 

j 


V 4 ax ~' 2 


dt 


(42) 


Furthermore, writing t = / max sin a, we get for the last two integrals: 


? 2 maxCOS 2 q 


I* -Xl-ra*-'_ dt = f_ 

J (i 2 -^ 2 -* 2 ), 2 ) 3 -U^-frVmaxSin 2 ^ 


da 


and 


‘ma: 

J 


v: 


,2 -r 2 

i max r 


^max COS 2 Ctf 


Q {b 2 +(a 2 -b 2 )rf [b 2 +(a 2 ~b 2 )t 2 max sin 2 af 


da 


Similarly, the second integral in equation (39) can be expressed as: 


(43) 


(44) 


V cos^-y/sin 2 ^ max -sin 2 ^ 
0 (* 2 +(a 2 -i 2 )sin 2 y9 ) 2 



7T 


V? 2 max~r £ _ j*_ f iLx c OS 2 Q- _ 

[b 2 +(a 2 -b 2 )t 2 J q (t> 2 +(a 2 -t» 2 )/ 2 max sin 2 a) 2 


da. 


(45) 


The last integrals described in the equations (43) (44) and (45) are evaluated using the 
formulas 2.563-1, 2.563-2, and 2.563-3 from [Ref. 5]. 

The marginal p.d.f. in azimuth is then given by: 


f<j> b {<t>b )~h( ( l ) b)~h {0b )’ sin (—o)-^- sin ! fe) 


where I (<t> )= —— ( l > bi a2 [b> 2 -a 2 ^b 2 [a 1 -D 2 sin 2 4) (3a 2 -4b 2 V 2 +(4fr 4 -2a 2 fr 2 V~a V 

3V rjr[a 2 -b 2 j 


and = (^ 2 ) , 

Z V r 0 / 'ly j — / 9 T > 

Vr(« -& j 


(46) 

(47) 

(48) 
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with y = a 2 sin 2 P max + b 2 cos 2 /3 max . 

E. MARGINAL AOA P.D.F IN AZIMUTH PLANE FOR THE SPECIAL 
CASE OF THE SEMI-MINOR AXIS OF THE SPHEROID GOING ZERO 

In this section we examine the special case of marginal AOA p.d.f., given by 

equation (46), for values of b^O. In order to evaluate the p.d.f. for this case the 

following simplifications are made on equations (47) and (48): 

y — a sin /? max + b cos /? max = b cos /? max (l *- an Amax) 

where tan/? max is given by equation (29). Since b^O the last equation can be 

L2 r\2 2 j. 

expressed as follows: y =- , . (49) 

D l -cr v J 

The terms y-b 2 and (y-b 2 ){p 2 -a 2 ) can be expressed as: 

y-b 2 ^b 2 ^^>0 a nd(y-b 2 )(D 2 -a 2 )=b 2 {a 2 -D 2 sm 2 ^ b ). (50) 


Plugging equations (47) and (48) to equation (46) the latter can be written as : 



(51) 

The term [(3a 2 -4 b 2 ^y 2 + [ 4 b 4 -2a 2 b 2 ^jy—a 2 b 4 jy can be expressed as: 

[(3a 2 -4b 2 )y 2 + fab 4 -2a 2 b 2 )y-a 2 b 4 ]j = 3a 2 y-2a 2 b 2 - a 2 ^--4b 2 y+ 4b 4 . (52) 
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The first three terms in the right side of equation (52) are second order terms, while the 
last two are fourth order terms with respect to b therefore can be neglected. Hence 
equation (52) reduces to: 

3 a 2 y - 2 a 2 b 2 -a 2 ^-4b 2 y + 4 b 4 = , aV . 

y (£> 2 -a 2 j 



(53) 

Combining equations (49), (50) and (53) and taking also in consideration which terms are 
second or fourth order with respect to ‘b‘, equation (51) can be written as follows: 




2 . pi 2 • 2 

■ (X + L) sin i 


+ D 2 cos 2 1 


^ 3D 1 cos 2 0 t 
p 2 -a 2 ) 


- 2 - P 2 -a 2 


D 2 cos 2 (i 


Carrying out some simplifications in equation (54) we get: 

h(h)~h(fa ) = ~D 2 sin 2 <j ) b ). 


(54) 


(55) 


Finally the marginal AOA p.d.f in azimuth plane of the spheroid for the values of b —» 0 
is given by: 

fhM = jj( a2 -D 2 sm 2 ^ b ) sin '(~-^)<<f> b < sin -1 (-g-). (56) 


F. DERIVATION OF MARGINAL AOA PDF IN AZIMUTH PLANE 

THROUGH A SECOND APPROACH 

As we have mentioned before the marginal AOA p.d.f in azimuth plane is given 

by the: 


Ami r h (fb’P) Amu r h(h-0) 

f<p b (0b ) = V J J r l sin Pd r bP = jr j jr b r b cos pdr b d(3 . 

“"Anax ~Anax r b\ {Qb'fi) 


(57) 
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Figure 9. The 2D coordinate system used for the derivation of the azimuth p.d.f. 

Figure 9, shows the intersected region between a plane at a specific (j) b value and 
the scatterer region volume. In the plane ^constant, we adopt a 2D Cartesian 
coordinate system x t yi. The terms of the integral in the equation (57), r b cos /? and 
r b dr b df3 can be written as: 

r b cos 15 = D cos <p + X cos a and r b dr b df3 = dx x dy { = dx 2 dy 2 . (58) 

In the x 2 y 2 2D Cartesian coordinate system as shown in Figure 9, the term /l cos a = x 2 . 
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According to the above equation (58) equation (57) can be written as: 


/?max r b 2 ^P) ^/a 2 -Z) 2 sin 2 ^ -Ja 2 -D 2 sm 2 0 b 

f<p b (0b)=v J | r b cos Pr b drdfi = A. J J(Z)cos^ + Acosa)dx 2 dy 2 

~dmax r b[ 0 0 


(59) 

The integral in equation (59) with respect to polar coordinate system A and a can be 
expressed as: 

2# Anax 

f<i> b (0b) = i J J (Dcos</> b + X cos a)MMa . (60) 

a=0 0 

Taking in consideration that the term A cos a is equal to X 2 , since X 2 is an odd function, 
equation (60) can be written as: 

271 Anax 

f</> t (0b) = v f J (A cos (p b + x 2 )AdAda = A jj DcosfoAdAda (61) 

a=0 0 

where \(0 b ) is the area of region R. 

Equation (61) is simplified as: 

A(&) = ££>cos^ \\MAda = ^-A 0b {(p b ) (62) 

where A t j >b (</> b ) = 7r-^(a 2 -D 2 sin 2 Substituting the above term in equation (62) we 

get for the marginal AOA p.d.f. in azimuth plane: 

fh^b) = ^^ L {a 2 -D 2 sm 2 <t> b ) sm~\-^)<<l> b <sm-%). (63) 

The last solution approach shows that the azimuth angular distribution is independent of 
b and coincides with equation (56). Hence, we conclude that the limiting equation (56), 
valid for b —> 0, is also true for all b. 
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Finally, using the closed form expression for the p.d.f. we evaluate the r.m.s 
angular spread in closed form: 


$rms ^0std 


where 


sin *(i) 

<Pld = = ( sin_1 te)F +f/J-i sin_1 fe)( 4 -^)-f+f^- 


sin 


~'h) 


(64) 


(65) 


G. MARGINAL AOA P.D.F. IN ELEVATION PLANE FOR THE SPHEROID 
SPATIAL CHANNEL MODEL 


As noted in section C of the present chapter, the marginal AOA p.d.f. in elevation 
plane is given by: 

«U£) 

ffi (P) = V J J r l cos M =W J k 3 {<Pb > P) - r l {<f>b > £)]cos Pd(t> b (66) 

MM *,{« 


0*2 (^h-^max 

fp(0)=w ji r bM’P)-^{(p b ,^)]cos J3d(/) b , \P\< p M (67) 

o 


The maximum elevation angle is obtained from the values that bound the 


spheroid in the xz-plane and can be found analytically from the equation (29) by setting 
the azimuth angle ^equals to zero: 


Pm =ta n _1 (^rj)- (68) 
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Using the equations (31) through (35), the integral in equation (67) can be simplified as 
follows: 


“max 

fp{P) = ir J cos ^k 


M 


+\ % 





rmax /— -.- / 

_ 2cos/? f y£~(/ Qcos 2 fl & -4A(;ff)C [B 2 (;0 )cos 2 ^ 

J *(£> \ A 2 (0) 

0 



^max 

_ 2g 2 (yg)co Sj g f 

34 3 (£)V J 

0 


cos 2 


</> b p 2 (/3)cos 2 <f> b -4A(j3)Cd</> b - 


2Ccos ft 
3A 2 (p)V 


“max 

J ■yjB 2 (fi)cos 2 0 b -4A(fi)Cd</> b 


(69) 


Substituting the functional forms of A(/3),B( ft) and C from equation (26) we get: 


■/>(£) = 


8Z>Vcos 3 


l{a 2 sin 2 0+b 1 cos* ffiv 


rrnax ,—- 

^- 3 - J cos 2 <j) b ^B 2 [P)cos 2 (j) b -4A{0)Cd(j) b - 


2 ( 2 2 \ 0max r—-- 

■ (70) 

0 

The upper limit for the azimuth angle ^ max , for a fixed ft, is the largest value of angle 
Qb that bound the region R , which is generated by the intersection of the cone 
P~ constant with the spheroid. In general, the intersection will include all the points 
r b 1 - r b — r b 2 ’ where and r b ^ are given in equation (27). Once again, at the boundary 
points it yields = r b ^. Thus 0 max is obtained by setting the discriminant 
) = B~(fi) cos 2 <j) b - 4A(j3)C equals to zero and solving with respect to (f> b : 

Anax = Sin- 1 tVl-(^) r } (71) 
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In order to be able to solve the integral in equation (70) we have to simplify the 
expression for the discriminant A(^) = B 2 (p) cos 2 <j) b -4A(J3)C under the square root. 
This is done applying the same approach as in previous section D. The discriminant, at a 
specific value of becomes equals to zero when the azimuth angle <p b coincides with 

the maximum value of <p b = <p max . Thus writing the discriminant as 

Afo) = Afo)— 

= [B 2 (p)cos 2 0 b - 4A(P)C] - [B 2 (P) cos 2 </> max - 4 A(j3)C ], becomes 

) = Wb A cos 2 P cos 2 <j> b - 4b 2 (j D 2 - a 2 \a 2 sin 2 /3+b 2 cos 2 /?)j- 

- [4 D 2 b 4 cos 2 P cos 2 <p m3x -4b 2 (p 2 - a 2 \a 2 sin 2 P + b 2 cos 2 p)\ (72) 

which gives: 

A {0b ) = A (0b ) - A (^max) = 4£> V cos 2 p(sin 2 0 max - sin 2 ) . (73) 

Substituting equation (73) into equation (70) gives: 


, t6 r ? f, 

P ?>(a~ sin 2 /3+b 2 cos 2 /3) l V 


'max i --- 

J cos 2 <p b (sin 2 <p mzx - sin 2 </> b )d<p b - 
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kind. Indeed with the help of the formulas 2.595-1 and 2.592-2 of [Ref. 5], these 
integrals take the form: 


r max 

sin^ J cos 2 <p b . 




max 


sin <p„ 


J 


(77) 


and 


sin ^ J $--zfc si * 2 h} ti h = £ ( sin $ max )~(sin^-l)K(sin^ max ) (78) 

max 0 

where ^(sin^^) and £ , (sin^ max ) denote the complete elliptical integral of the first and 
second kind respectively [Ref. 5]. 

Substituting the results from the equations (77) and (78) into equation (74) we get 
the final expression for the marginal AOA p.d.f., in elevation fp(fi) in a closed form as 
follows: 


6(tf 2 sin 2 fi+b 2 cos 2 fifv [(^ + S ^ n ^max )£(sin 0 max ) COS ^ max AT(sin ^ max ))] 

~ ^ (sin ) - cos 2 0 max ^(sin 0 max )], \p\ < P M (79) 


with P M = tan Sin^ = ^l-(^) 2 ~ and V =jm 2 b. 
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III. NUMERICAL RESULTS AND DISCUSSION 

Having derived closed form expressions for the marginal AOA p.d.f. in azimuth 

and elevation planes, we now present in this chapter some numerical results for each 
marginal AOA p.d.f.. Recall that D is the distance between transmitter and receiver, a is 
the semi-major axis of the spheroid, and b is the semi-minor-axis of the spheroid. 

A. PLOTS FOR MARGINAL AOA PDF IN AZIMUTH PLANE 

In this section we consider the cases of ‘D’, ’a’ and ‘b’ as given in Table 2 and 
generate some plots: 


D(m) 

a (m) 

b(m) 

Corresponding 




Figure 

1000 

100 

0.001 

10 

1000 

100 

50 

11 

1000 

200 

150 

12 


Table 2. Parameters of D, a and b used in plots for marginal AOA p.d.f. in the azimuth 
plane. 

Figure 10 has been plotted in order to check how the azimuth angular distributions 
of the spheroid and the circular scattering models are related to each other. The circular 
scattering model is presented in [Ref. 3]. The angular distribution for the circular model 
is given by equation (80): 
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(80) 


f<p b (0b)~ 


2D cos 


D 2 cos 2 c 


> -D 2 +a 1 


m 




0 ,else 



Figure 10. Plot of (<p b ) versus^ forD=1000m, a=100m and b=0.001m. Also shows 
the AOA p.d.f. for a circular spatial model for D=1000m and a=100m. 


In Figure 10, we plot the p.d.f. for a/D= 0.1 and b/D= 10' 6 . Comparing the plot of the 
spheroid model for small values of ‘b‘ with this of the circular model, we can observe that 
there exists a significant difference between them. 

In Figures 11 and 12 we plot the p.d.f. for a/D= 0.1, 6=50m and for a/D= 0.2, 
b- 150m using the two different formulas for the angular distribution in the azimuth 
plane. 
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Angle of Arrival in Azimuth Plane(degrees),<}> b 


Figure 11. Plot of (<j> b ) versus <j) b for D=1000m, a=100m and b=50m 
(using the two different formulas). 



Angle of Arrival in Azimuth P!ane(degrees),<j> b 


Figure 12. Plot of (<j> b ) versus <p b for D=1000m, a=200m and b=150m 
(using the two different formulas). 































As shown Figure 11 and 12, the two different formulas give identical results for the 
angular distribution in azimuth plane. Additionally, we conclude that the higher the ratio 
a/D, the lower the probability for the received signals at the base station to be restricted in 
a small angular region centered about the axis which connects the base and the mobile 
stations. 

Finally in Figure 13 we plot the r.m.s angle spread in azimuth plane as given in 
equation (64) versus the ratio a/D. 



Figure 13. Plot of the r.m.s angular spread versus a/D. 


Figure 13 shows that for values up to 0.5, of the ratio a/D, the profile of the r.m.s angle 
spread is almost linear. Note that the r.m.s angle spread in azimuth plane is independent 
of b as proved with the derivation of the second formula for the p.d.f in azimuth plane. 
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B. PLOTS FOR MARGINAL AOA PDF IN ELEVATION PLANE 


D(m) 

a(m) 

b(m) 

Corresponding 

Figure 

1000 

100 

1 

14 

1000 

100 

30 

15 

1000 

100 

50 

16 

1000 

200 

50 

17 

1000 

500 

50 

18 


Table 3. Parameters D, a, and b used in plots for marginal AOA p.d.f. in the 
elevation plane. 


Figure 14 is generated to demonstrate the angular distribution in elevation plane 
for small values of b. 
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Angle of Arrival in Elevation Plane(degrees),{3 

Figure 14. Plot of fp(fi )versus for D=1000m, a= 100m 

and b=lm.. 

As we can see form the Figure 14 the angular distribution in elevation plane for small 
values of b looks like a dirac function as it is expected. 

Figure 15 and 16 are plotted for demonstration purposes, giving information of 
how the changes of the ratio b/D for a given ’a’ affect the angular distribution in elevation 
plane. 
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Angle of Arrival in Elevation Plane(degrees),p 

Figure 15. Plot of fp(fi) versus (5 for D=1000m, a= 100m 
and b=30m.. 

Figure 15 shows the the angular distribution in elevation plane for a ratio b/D= 0.03 and 
a/D= 0.1. 



Angle of Arrival in Elevation Plane(degrees).f} 


Figure 16. Plot of fpifi) versus for D=1000m, a=100m 
and b=50m. 
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Figure 16 shows the p.d.f for a ratio b/D= 0.05 and a/D=0.l As seen, compared to Figures 
15 and 16, the higher the ratio b/D the lower the probability to receive multipath 
components from a small region fi BW at the base station. 

Figure 17 is generated, as set with Figure 16, in order to show the effect of the 
length of the semi- major axis a on the angular distribution. 



Angle of Arrival in Elevation Plane(degrees),p 


Figure 17. Plot of fp{fi) versus /? for D=1000m, a=200m 
and b=50in meters. 

Figure 17 shows the p.d.f for a ratio b/D= 0.05 and a/D— 0.2. Comparing Figuresl6 and 
17, it is seen that p.d.f in elevation plane is almost independent of the length of semi¬ 
major axis of the spheroid for small values of the ratio a/D. Figure 18, is generated to 
show the effect of higher values of the a/D ratio on the angular distribution. 
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Angle of Arrival in Elevation PIane(degrees),p 


Figure 18. Plot of fpifi) versus forD=1000m, a=500m, 
and b=50m. 

Comparing Figure 17 and 18, it was observed that there exists a difference between the 
two angular distributions. 

Finally in Figure 19 we plot the r.m.s angular spread of the angle of arrival in 
elevation plane, as given by the equation (81), versus the ratio b/D for specific values of 
a/D. 



(81) 


The mean value of the angular distribution at the base station is zero, and the integral is 
evaluated numerically. 
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Figure 19. Plot of the r.m.s angular (3 rms spread versus b/D 
for various ratios a/D,for the spheroid scattering model. 


As shown in Figure 19, there exists a linear relation between the r.m.s angular spread and 
the ratio b/D independently of the ratio a/D. 

Eventually for validation purposes of the two derived, in this thesis, closed 
forms of the angular distributions, we plot in Figure 20 both p.d.f. with a almost equals to 
b (b is not made exactly equal to a to avoid dealing with limiting process). 


42 

















$ b ,p in (degrees) 

Figure 20. Plot of versus <t> b and fp{0) versus ft for D=1000m, 

a=100m and b=99.99m and for D= 1000m, a=100m and b=100m respectively. 

Figure 20, shows the angular distribution in azimuth and elevation plane for 
a/D=b/D-0. 1. From the plot we can see identical results for both angular distributions, as 
expected due to the symmetry of the model. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

In this thesis we proposed a geometrically based single bounce scattering 3D 
model. It has been assumed that the scatterers are distributed about the mobile with a 
constant density function in a volume as defined by the geometry of the spheroid. The 
distance between the base station and the mobile station is D. The focus of the effort of 
this research was in the derivation of the mathematical expressions which described the 
spatial statistics properties of the received signals in the presence of the three-dimensional 
multipath scattering. 

After presenting a general approach for the derivation of the joint TOA/AOA and 
marginal AOA p.d.f.s in terms of an arbitary scatterer probability density function, we 
presented the derivation of the angular distribution in azimuth and elevation plane for the 
spheroidal model in closed expression forms.. 

Making use of these quantities we studied the behavior of the angular spread for 
each p.d.f as function either of the spheroid dimensions or of the distance between the 
transmitter and the receiver. For either angular distribution the model results in a higher 
probability of arrival along the line of sight direction. Also, we observed that the higher 
the ratio of each spheroid’s dimension over the distance D the lower the probability that 
the multipath components are confined to a small region. An excellent agreement exists 
in the plots which are generated from the two different formulas of the p.d.f in azimuth 
plane. 
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Analytical and numerical calculations were performed for the r.m.s angular spread 
in azimuth and elevation planes respectively. The effect of the ratio a/D or b/D on the 
corresponding angular spread was investigated. 

B. RECOMMENDATIONS 

Further research is required to enhance this model. It will be a interesting to derive 
the marginal time of arrival p.d.f or to investigate the effect of the angular spread on the 
correlation observed between the received signals an antenna array. Another area of 
interest should be the validation of this model. This could be done comparing with 
measurements taken over similar macrocell environments. 
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